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We consider theoretically an armchair Carbon nanotube (CNT) in the presence of an electric 
field and in contact with an s-wave superconductor. We show that the proximity effect opens up 
superconducting gaps in the CNT of different strengths for the exterior and interior branches of 
the two Dirac points. For strong proximity induced superconductivity the interior gap can be of 
the p-wave type, while the exterior gap can be tuned by the electric field to be of the s-wave type. 
Such a setup supports a single Majorana bound state at each end of the CNT. In the case of a weak 
proximity induced superconductivity, the gaps in both branches are of the p-wave type. However, 
the temperature can be chosen in such a way that the smallest gap is effectively closed. Using 
renormalization group techniques we show that the Majorana bound states exist even after taking 
into account electron-electron interactions. 

PACS numbers: 73.63.Fg, 74.45. +c 



Introduction. Majorana fermions in solid state s 
terns have attracted considerable attention recently 
Ho| . In particular, the possibility of realizing them as 
bound states at the ends of semiconducting nanowires in 
the proximity of an s-wave bulk superconductor has led 
to much activity. Such setups require a Zeeman split- 
ting, typically generated by an external magnetic field 
, that must be larger than the proximity induced gap 
to induce an effective p-wave superconductor in the topo- 
logical phase. Such a magnetic field, however, tends to 
destroy the gap in the bulk superconductor itself, and 
thus a delicate balance must be found [12| . It is therefore 
very desirable to search for Majorana-scenarios which do 
not require magnetic fields. 

One of the prerequisites for a Majorana bound end 
state (MBS) is the existence of helical modes, i.e. modes 
which carry opposite spins in opposite directions. It has 
been shown recently that such helical states are induced 
in Carbon nanotubes (CNT) via spin-orbit interaction 
(SOI) by an external electric field 13|, [14J. This mecha- 
nism works optimally for a special class of metallic CNTs: 
armchair CNTs (N, N). This class is characterized by a 
spin-degenerate low-energy spectrum around the two in- 
equivalent Dirac points, K and K' . This degeneracy can 
be lifted by an electric field which gives then rise to he- 
lical modes. 

However, when putting the CNT in contact with an 
s-wave superconductor (see Fig. [IJ with the goal to gen- 
erate MBS the following problem is encountered. The 
two Dirac points K and K' are Kramers partners (see 
Fig. [2]) and thus the superconducting pairing induced via 
the proximity effect will involve both of them, i.e. left 
(right)-moving electrons from the branch at K get paired 
with the right (left)-moving electrons from the branch at 
K' to form an s-wave Cooper pair with zero total momen- 
tum. This then immediately implies that there will be 
two superconducting gaps, an 'exterior' one, A e , and an 
'interior' one, A;. Thus, in general, we expect two MBS 
at each end of the CNT (i.e. four in total). This, how- 




FIG. 1. (a) An armchair nanotube (cylinder) is placed on top 
of a superconductor (blue slab). The x-axis points along the 
nanotube. An electric field E is applied perpendicular to the 
nanotube in the y-direction. There are two non-equivalent 
lattice sites: A (light red) and B (light green), (b) The dis- 
tances between the superconductor surface and the atoms of 
sublattice A (dark-red row) and of sublattice B (dark-green 
row) are assumed to be the same. Thus, the tunneling ampli- 
tudes to the different sublattices are (nearly) equal. 



ever, is problematic as the Majorana pair at a given end 
can easily combine to form a single fermion by some local 
perturbation. Thus, the question then arises whether it 
is possible to find a regime with only one MBS at each 
end 15j. As we will show, the answer is affirmative but 
under rather stringent conditions. One of them requires 
a comparable tunnel coupling of the A and B sublattices 
of the CNT to the superconductor, see Fig. Q] Using the 
interference mechanism first described by Lc Hur et al. 
[lfjj ]. we will show that for this particular case the exte- 
rior (interior) gap A e (Aj) gets enhanced (suppressed) 
due to constructive (destructive) interference in the tun- 
neling process. If A e u, is smaller (larger) than the gap 
opened by the electric field, then the coupling between 
the two Dirac points is of p-wave type (s-wave type) . This 
leads to two regimes for MBS. The first one involves a 
scenario wherein only one of two branches has a p-wave 
gap, thus giving rise to only one MBS at each end of the 
CNT. In a second regime, where both branches have a 
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potential for p-wave pairing, the temperature T can be 
choosen to lie between the two gaps A e and Aj, so that 
only the exterior branches will go fully superconducting, 
whereas the interior branches stay normal. Again, a sin- 
gle pair of MBS in the CNT emerges. We further in- 
vestigate the effect of interactions on the MBS. This is 
particularly important for the second regime due to the 
presence of gapless states from the interior branches that 
could be harmful to the MBS. However, using bosoniza- 
tion techniques we will dispel these concerns and show 
that for screened interactions the MBS remain stable al- 
though they can get substantially dclocalized similar to 
the simpler case of Rashba wires || . 

Nanotube spectrum. We consider an armchair CNT 
in the presence of an electric field E applied perpen- 
dicular to the CNT axis (see Fig. [T]). Taking into ac- 
count the spin-orbit interaction the low-energy sector is 
described by an effective Hamiltonian around the Dirac 
points given by [IH 



H = hvFkT 3 a 2 + T 3 eE£S z <r 2 + aS x ai, 



(1) 



where k is the momentum along the nanotube axis taken 
from the Dirac point, oi is the Pauli matrix on the sub- 
lattice space (^4, B) associated with the honeycomb unit 
cell, and Si is the spin operator with eigenvalues ±1. 
The Pauli matrix Tj acts on the K, if'-subspace. Here, 
vf — 10 6 m/s is the Fermi velocity, and the parameter 
a arises from the interplay between SOI and curvature 
effects ID, 17 , 18| . In the framework of the tight-binding 
model, a = — 0.08 meV/i?[nm], where R is the radius of 
the CNT [13]. The parameter £ ~ 2 x 10" 5 nm is given by 
a combination of hopping matrix elements, on-site dipolc 
moment, and SOI [141 ]. 

The spectrum given by H (Eq. []} consists of four 
branches (see Fig. [J) 



e n (k) = ±eE£ ± yV + (hv F k) 2 



(2) 



for each Dirac point. In the following, we label the four 
branches by n = 1, ...,4. For each k, n = 1 corresponds 
to the highest eigenvalue and n = 4 to the lowest. The 
remarkable feature of the spectrum is the existence of he- 
lical modes, which carry opposite spins in opposite direc- 
tions. The average value of the spin in the direction par- 
allel to the axis of the CNT {(S x )) or parallel to the ap- 
plied electric field {{S y )) is equal to zero. The projection 
of the spin along the z-direction is equal to (>!5 2 ) = sin£, 
where £ is defined by £ — &ics,uv(hv pk / \J a 2 + (hvFk) 2 ) 
and depends on the wavevectror k. Note that the eigen- 
vectors ip e Jx an d ^Jk' are independent of the electric 
field E. For a (10,10)-CNT in an electric field of strength 
lV/nm, and with a Fermi level /i tuned between the two 
lowest electronic states polarizations close to 90% can be 
reached fl4j |. 

Proximity induced superconductivity. If a CNT is in 
contact with an s-wavc bulk superconductor, then the 
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FIG. 2. The energy spectrum around the Dirac points K, K' 
for a (10,10)-CNT in an electric field E = 1 V/nm, which con- 
sists of exterior (full line) and interior (dashed line) branches. 
Each branch of the spectrum is characterized by the sign of 
the spin projection along the z-axis {S z ) (red: spin down, 
blue: spin up). The Fermi level fi lies inside the gap given by 
2eS£, and S = e££ + a - u. 



proximity effect induces superconductivity also in the 
CNT which at the BCS mean-field level is described by 



E 



(A dC ; 



f ct 
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(3) 



where we concentrate on the contribution coming from 
the 7r-bands formed by the radial p r -orbitals 14, [l9| . 



Here, 



„(t) 

"ip r s 



are the standard fermionic annihilation (cre- 
ation) operators, with s and s denoting opposite spin 
states [14{ . The sum runs over atoms which are in con- 
tact with the bulk superconductor: i and j belong to 
the same sublattice, whereas i' and j' belong to different 
sublatticcs. Gcnerically, one can assume that the lattice 
constant of the superconducting material is not commen- 
surate with the one of graphenc. The CNT is placed in 
such a way that the distance from the superconducting 
surface to the A and B atoms is the same (see Fig. 1), 
which is satisfied for the special case of armchair CNTs. 
This ensures equal probability amplitude for tunneling to 
either sublattice. Since the phase of the superconducting 
order parameter A^/„ can be chosen arbitrary, we assume 
them to be real. The coupling terms in Eq. ^ conserve 
momentum, so they pair Kramers partners from the op- 
posite Dirac cones. The process in which the Cooper pair 
tunnels from the superconductor to either one of the sub- 
lattice a is written as 



(4) 



where the indices k and R denote opposite Dirac points. 
The operators ifi asK and Ci PrS are connected via Fourier 
transformation [14| . The pairing term between electrons 
in different sublattices are 



iA n ^2 sgn(s)ipl SK tljl gR + h 



(5) 
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To simplify the notation we introduce Pauli ma- 
trices rji which act on the particle-hole subspace, 
and we work in the basis (\&, with \& = 

This allows us to rewrite Eqs. Q and (|5|) in a compact 
form 

H sc = —T)2TlS z Ad <7 3 + T]lTlS z A n £71. (6) 

In the same particle-hole basis, the effective Hamiltonian 
Eq. ([!} can be rewritten as 

H = 7}{frvFk(T2 + im T 3S z (J2 + ar) S S x ai). (7) 

To express the coupling between the different energy 
states in a canonical form we work in the basis of eigen- 
vectors {^nK^nK^nK'^liK'}- For tnG states at the 
Fermi level, n = 2, H sc takes the form 

M^2K4 l 2K-^2K^K') + h.C, (8) 

I— e,i 

with different coupling strengths for the exterior (e) and 
interior (i) branches, 

A e/i = A d ±A n |sinC|. (9) 

We note that the sign reflects the constructive and de- 
structive interference, resp., in the tunne ling process from 
the bulk-superconductor into the CNT [16(. We neglect 
the term A n cos ( characterizing the coupling between 

ip2K an d V'Ik' which are separated by the particle-hole 
gap 2a, see Fig. [2] In the following we consider the 
limit of equal diagonal and non-diagonal parameters, i.e., 
Ad = A n . We note that for k ^> a/HvF the coupling 
between the interior branches is close to zero and that 
between the exterior branches is equal to 2 Ad- We show 
that this asymmetry in the coupling strengths is crucial 
for the existence of Majorana bound states in CNTs. 

Majorana bound states. Next, we obtain the MBS fol- 
lowing the derivation of Ref. |8|. For illustrative pur- 
poses we derive the bound states that arise by con- 
sidering the exterior branches. The field correspond- 
ing to the exterior branch is defined as, tp e (x) = 
ip2\(x)e t( - kF+K ^ + ip% K ,{x)e~ l( * kF+K ^ x , where ^f K {x) 
and ip2K'( x ) are the slowly moving right and left com- 
ponents about the K and K' points, resp. Denoting 
the length of the CNT by L, the boundary conditions, 
ip e (% = 0) = i/) e (x = L) = 0, yield the restriction 
^2k( x ) = ~ ^2K'(~ X )- Thus, the kinetic term is given 

by Hq 1 ^ = —wf J-T l i ) 2K{ x )^^2k{ x )^ an d the p-wave 
pairing term between the exterior branches by 

- A e J dxsgn(x)[^ K (x)^ K (-x) + h.c] . (10) 

Solving for the zero energy mode localized around x = 0, 
we obtain the MBS 

*f (x) (xj e sm[(K + k F )x]e- x/ ^, (11) 
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FIG. 3. The particle-hole spectrum of a CNT (10,10) in the 
presence of an electric field E with the Fermi level tuned inside 
the energy gap between the two upper branches (solid black 
lines). All energies are counted from the Fermi level fi = 
0.11 meV (see Fig. [2]). By proximity effect, superconducting 
gaps A e ,i are opened at the Fermi points, (a) We keep E 
constant at 1 V/nm and change the strength of the induced 
superconductivity. For Ad = 5 fj,eV < A c i both exterior and 
interior branches are in the p-wave phase (dotted blue line). 
At the critical value A^ = 23 fieV = A c i the gaps induced by 
the proximity effect and by the electric field are equal (dot- 
dashed red line). For Ad = 30 /j,eV > A c i only the interior 
branch is in the p-wave phase (dashed green line). Keeping 
Ad constant at 11 \ie~V and changing E, one goes from a 
regime [dashed blue line in (b)] with E — 0.4 V/nm < E c i = 
0.6 V/nm in which only the interior branch is in the p-wave 
phase to a regime [dashed green line (c)] with E = 1 V/nm > 
Ed = 0.6 V/nm in which both the exterior and the interior 
branches are in the p-wave phase. 

where j e = 7J, and it is assumed that the localization 
length, given by £ e = Sui?/2A e , satisfies £ e -C L. Simi- 
larly for the interior branches, with the index e replaced 
by i. 

In general, the Majorana modes arising from the inte- 
rior and exterior branches at the same end of the CNT 
are not protected and can combine into a fermion. To 
avoid such a scenario one needs to ensure the presence of 
only one single MBS at each end of the CNT. This can 
be achieved in two ways. 

First, there is a window where the electric field E can 
be chosen in such a way that the superconductivity in the 
exterior branch can be tuned into the s-wave phase, while 
the interior one still remains in the p-wave phase. In this 
case, only the interior branch supports a MBS at each 
end of the nanotube, and we refer to this as a topologi- 
cal phase of the CNT (see dashed green line in Fig. [3^,). 



4 



Concretely, such a regime is reached for A e (kp) > 8 > 
Aj(/cf), where kp « ^/(/i + eE£) 2 — a 2 /hvp and 6 = 
eE£ + a — (i. With Eq. (j9|) this criterion becomes equiv- 
alent to A c2 > Ad > A c i, where A cl / c2 = 6/(1 ± sin £). 
For a given value of Ad, the experimentally viable ap- 
proach to drive the system into the topological phase is 
to tune the electric field E. Indeed, for E c \ > E > E c2 
(see Fig. [3Jd) the exterior branch is in the s-wave phase, 
while the interior one is in the p-wave phase. The crit- 
ical value of the electric field E c i (E c2 ) is determined 
by the condition S = A e (kp) (S = A,-(fcp)). Similarly, 
we can tunc the Fermi level. [In passing wc note that 
the gap eE£, and thus 6, get enhanced by interaction ef- 
fects around k = [1^ , which is useful for experimental 
realizations. However, for simplicity wc will ignore this 
feature here.] 

Second, in the regime A^ < A c i (see dotted blue line in 
Fig-Eb) or E > E cl (see Fig. [3b ) both branches are domi- 
nated by p-wavc pairing. If the temperature is lower than 
both gaps, i.e. ksT < A e .j, then there is an even number 
of MBS at each end of the nanotube, and the CNT is in 
the topologically trivial phase. However, in the interme- 
diate regime with A e > ksT > Aj, the interior gap Aj 
is closed and the Majorana states are removed, yet those 
from the exterior branches remain, and the CNT is again 
in the topological phase. In the following we consider 
this latter scenario and discuss the role of interactions 
coming from the gapless states of the interior branch. 

Interactions. Interactions effects arc most conveniently 
described by linearizing the spectrum of the fermionic 



fields ip2K an d near the Fermi momentum and ex- 
pressing them in terms of the bosonized fields. The 
quadratic part of the bosonized Hamiltonian thus ob- 
tained has the following form, 

n ° = \ E {v n K n (d x 6 n ) 2 + ^(&0„) 2 }, (12) 

Z n=± A " 

where d x (j)+ and d x 4>- are the sum and difference of den- 
sities between the two fermionic bands. The fields conju- 
gate to them are defined as, + and resp. The param- 
eters K+ ~ 1-Uo/irvaxidK- ~ l + {l-{S z ) 2 )U 2kF /^v 
encode information about the interactions and the renor- 
malized velocities are given as v + ~ Vp+Uojjr, and w_ ~ 
v F + b'(l + (S z ) 2 )/Att, where the b'-term [Hi, HI] is due 
to the backscattering contribution. Here, Uo,2k F denotes 
the Fourier component of the screened Coulomb interac- 
tion. Since (S z ) 2 < 1 and thus AT_ > 1, wc conclude [23| 
that the forward scattering term oc J dxdr cos(v / 87T0-) 
scales to zero. 

Additional terms induced by the proximity effect lead 
to a modified Hamiltonian given by 



U = •Ho + TT-^COS 

Zira 



2tt(6L| 



•>-)+- — COS 



<2tz{ 



Since we assume here A e > ksT > Aj, the second term 
due to the interior branches ( A,) is smeared out by tem- 
perature effects and will not be considered. 

Using standard techniques 0, [24j , we derive the fol- 
lowing renormalization group (RG) equations, 

4 
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(14) 



X_A%(l + 7 ) 2 
7(1 - l)K+ 



4 (\+i){K+K- + l) ' 
1 K- ~ 



f 2 



2K 4 



(15) 
(16) 
(17) 



(13) 



where the flow parameter I = ln[a/ao], / = 2A e a, and 7 
is the ratio of the velocities v + jv_ . We note that for the 
non-interacting case 7 is already at its fixed point, 7=1, 
and including interactions (the repulsive interactions are 
assumed to be well screened) causes only a small devia- 
tion from unity 25, 2(|. Thus, it is convenient to assume 
7 = 1, and under this assumption A + A_ is a constant, 
given in leading order by unity. Above RG equations 
now acquire the simple form dR/dl = / 2 /2 and df /dl = 
/(2-1/iJ), where R= (l/A + + A_)/2. These equations 
are exactly the same as in Ref. Q derived for interact- 
ing spinless fermions in an effective p-wave regime. We 
conclude that for a CNT whose initial values of the pa- 
rameters lie in the regime fa > 2\j2R^ — ln(2i?oe) has 
its RG flow such that both K + and K- approach the 
non-interacting value. At this point the problem can be 
refcrmionized into a set of decoupled gapped and gapless 
fermions and for a strongly screened CNT with initial 
value e.g. K + = 0.8 the localization length £ e increases 
by 25%. Therefore, wc conclude that the MBS which 
arise from gapped fermions remain protected even in the 
presence of interacting gapless fermions and simply ac- 
quire a rcnormalizcd £ e . 

Conclusions. We have shown that an armchair CNT 
with helical modes generated by an external electric field 
is a promising candidate material for Majorana bound 
states. By placing the CNT on top of an s-wave super- 
conductor and tuning the Fermi level and the electric 
field, one can induce pairing of Kramers partners from 
opposite Dirac points. This pairing opens up inequiva- 
lcnt gaps for the exterior and the interior branches. The 
Majorana modes obtained are stabilized by cither tun- 
ing the electric field such that the exterior gap acquires a 
predominantly s-wave character or by increasing the tem- 
perature to remove the pairing in the interior branches. 
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